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A RAY THEORY FOR NONLINEAR SHIP WAVES 
AND WAVE RESISTANCE 


Bohyun Yim 
David W. Taylor Naval Ship Research and Development Center 
Bethesda, Maryland 20084, U.S.A. 


Abstract 

Analytical and numerical methods for 
application of ray theory in computing ship 
waves are investigated. The potentially 
important role of ray theory in analyses of 
nonlinear waves and wave resistance is demon- 
strated. The reflection of ship waves from 
the hull boundary is analyzed here for the 
first time. 

When a wave crest touches the ship sur- 
face, the ray exactly follows the ship surface. 
When the wave crest is nearly perpendicular to 
the ship surface the ray is reflected many 
times as it propagates along the ship surface. 
Many rays of refiected elementary waves 
intersect each other. The envelope to the first 
reflected rays forms a line like a shock 
front which borders the area of large waves 
or breaking waves near the ship. 

Fort the Wigley hull, ray paths, wave 
phases, and directions of elementary waves 
are computed by the ray theory and a method of 
computing wave resistance is developed. The 
wave phase is compared with that of linear 
theory as a function of ship-beam length 
ratio to identify the advancement of the bow 
wave phase which influences the design of bow 
bulbs. The wave resistance of the Wigley hull 
is computed using the amplitude function from 
Michell's thin ship theory and compared with 
values of Michell's wave resistance. The total 
wave resistance has the phase of hump and hollow 
shifted considerably. 


Introduction 


Significant developments in the ship wave 
theory have been made in recent years. These 
include the application of ray theory!»2* and 
the experimental discovery of a phenomenon 
called a free surface shock wave. 


Because a ship is not thin enough to apply 
both the thin ship theory and the complicated 
free-surface effect, theoretical development 
of an accurate ship wave theory has been slow. 
The problem should be evaluated differently 
from the conventional means. Ray theory has 
been used in geometrical optics or for waves 
having small wave lengths. Ursell* used the 


*A complete listing of references is given 
on page 15. 


theory to consider wave propagation in non- 
uniform flow, and Shen, Meyer, and Keller? used 
it to investigate water waves in channels and 
around islands. Recently Keller“ developed a 
ray theory for ship waves and pointed out that 
the theory could supply useful information 
about the waves of thick ships at relatively 
slow speeds. However, he had difficulty ob- 
taining the excitation function for wave ampli- 
tude and solved only a thin-ship problem. Inui 
and Kajitani! applied the Ursell ray theory* to 
ship waves, using the amplitude function from 
linear theory. 


Because the ray is the path of wave energy, 
it is not supposed to penetrate os ship surface; 
and this was emphasized by Keller. However, 
neither Keller nor Inui and Kajitani consider 
non penernaragn of the ray seriously. Recently 
Yim° found the existence of rays which emanated 
from the ship bow and reflected from the ship 
surface. 


In the present paper, further study has 
revealed many bow rays reflecting from the 
ship surface, creating an area in which these 
rays intersect each other. The envelope to 
the first reflected rays forms a line like a 
shock front which borders the area of large 
waves or breaking waves near the ship. This 
will be referred to as the second caustic. This 
phenomenon was observed by Inui et al.’ in the 
towing tank and called a free surface shock 
wave (FSSW). Much has been done concerning the 
theoretical investigation of FSSW by introduc- 
ing a fictitious depth for a shallow water non- 
dispersive wave. “ 


The ray equation and the ship boundary 
condition are analyzed further to show the 
existence of reflecting rays, except in the 
case of a flat wedgelike ship surface. The 
ray path is very sensitive to the initial bound- 
ary condition near the bow or the stern, and 
should be identified at downstream infinity 
not by the initial condition. The conservation 
law of wave energy in the nonuniform flow is 
different from that of the uniform flow due to 
the exchange of energy with the local flow.8 
Therefore, the linear wave amplitude as a 
function of initial wave angle near the bow 
or stern is meaningless and cannot be used in 
the ray theory. It is shown that the amplitude 
function for the ray theory matched with the 


linear theory far downstream is reasonable and 
likely to produce a reasonable result as in 
the two dimensional theory. 


The wave phase in the ray theory is ob- 
tained, together with the ray path, indepen- 
dent from the amplitude function. Both the 
phase and the ray path are quite different from 
the prediction of linear theory near the ship. 
The rays far downstream are straight as in the 
linear theory yet have phases different from 
the linear theory; they are parallel to the 
linear ray with the same wave angle but do 
not coincide. The difference in wave phases 
is the main factor that makes the wave resis- 
tance different from the linear theory. Ray 
paths and phase difference are computed for 
various parameters of the Wigley ship, with and 
without a bulbous bow. The ray paths for 
different drafts and different beam-length 
ratio of the Wigley hull are slightly different, 
widening the wave area near the hul! for the 
wider beam and/or larger draft. However, the 
phase difference is more sensitive to the beam- 
length ratio and/or draft-length ratio, by 
always advancing the linear wave phase. 


The wave resistance of the Wigley hutl is 
computed and shown to have a considerable shift 
of phase of hump and hollow. 


The most interesting phenomenon of a ship 
with a bulbous bow is the reduction of slopes 
of rays and the second caustic, i.e., the 
larger the bulb, the greater the reduction. 
This fact was observed in the towing tank. 
There exists a bulb size yack totally elimi- 
nates the reflecting rays.! However, the phase 
difference due to the bulb!! is very small 
showing that the phase difference, which has 
been observed in the towing tank, is the effect 
of nonuniform flow caused by the main hull. 


Ray Equations 


The concept of ray theory in ship waves is 
analogous to the concepts used in geometrical 
optics and in geometric acoustics. A ship is 
considered advancing with a constant velocity 
-U which is the direction of the negative x 
axis of a right handed rectangular coordinate 
system O-xyz with the origin O at the ship bow 
on the mean free surface, z=0, z is positive 
upwards. 


First the phase function s(x,y,z) is de- 
fined so that the equation (s = constant) re- 
represents the wave front where the value of s 
is the optical distance from the wave source, 
e.g., the ship bow. When Keller? developed 
his ray theory of ships by expanding boundary 
conditions and a solution, which should satisfy 
both the Laplace equation and the boundary 
condietons) in a series of Froude number squares 
F*, he obtained: 


2 
WH 10 (1) 


2 
a. = -i (S, +Vo * Vs) at z=0 


By eliminating Sz from these two equations he 
obtained a dispersion relation 


2 2.35 2 e 
(s, + a, ) (Ss, + oy 5, + thy a) at Zz 0 
(2) 

where » is the double model potential. When 
steady state motion is assumed with respect to 


the moving coordinate system, 0 - xyz, 


Gee 0 


When the angle between the normal n to the 
phase curve s and the axis is denoted by 6, 


n = icosé + jsiné (3) 
Then the wave number vector is defined by 


Be Tefal) S ae a 


1 x y 
Pa = = (4) 
= nk = ik cos6 + jk sinO at z = 0 
From equations (2) and (4) 
1 2 
Us -(; cosé + v =| (5) 


where 


5oR =u and ats, =v 


These results are an approximation within 
the order of F4, and the phase function and its 
related equations are all limited to their 
values at z= 0. Thus, from now on, unless 
otherwise mentioned, all the physical values 
are at z= 0. The ray equation of ship waves 
is obtained from the irrotationality of the 
wave number vector 


ok dk 
ez ela (6) 


ax dy 
From Equations (4) through (6), 
{2 siné (u sin® - v cosé) 


+ cos8 (u cos® + v siné)} - 


+ {2 cosé (-u sin@ + v cosé) 


+ sin@ (u cosé + v sin@)} a Q) 
= 2 sind (cosa 22 + sind 2%) 
8 98 ax ax 
du ov 
- 2 cos@ (cosé Oa + sind —) 
y ay 


This is the ray equation which can be 
solved by the method of characteristics, and 
is equivalent to simultaneous ordinary differ- 
ential equations. 


dy * {v- 4 sin@ (u cosé + v siné)} (8) 


dx  {u- 5 cos6 (u cos® + v siné)} 


aG du av 
ae CalZ s6 — ahs 
qe ={2 sind (cose 5 + sind x) 


- 2 cos@ (cos8 of + siné Pw Mi2 siné (u sind (9) 


- v cos8) + cos@ (u cos + v siné) } 


Here u cos6 + v sin 6 is the velocity component 
normal to the phase curve of the flow relative 
to the ship. Because the wave is stationary 
relative to the ship, the phase velocity 
through the water surface should be 


- u cos@ - v sin#8 


The group velocity is one-half of the phase 
velocity. Thus the ray direction in Equation 
(8) is along the resultant of the group velo- 
city taken normal to the phase curve and of the 
velocity of the basic flow as Ursell has shown. 
Because the wave energy is propagated at the 
group velocity, the ray path is interpreted as 
the path of energy relative to the ship. This 
can be obtained by solving Equations (8) and 
(9) with the proper initial condition. The 
phase s can be obtained from equations (4) by 


e+ [kar 


as in potential theory; s(x,y) is a function of 
(x,y) but is unrelated to the integration path. 
However, 


ds = k cos@ dx + k sind dy (10) 


can be solved together with the ray equations 
along the ray path. 


Rays of Ship Waves and Linear Theory 


To investigate the path of a ray of a 
ship wave, a ship, represented by a double 
model source distribution m (x,y) on 
y = 0, h>z>-h, is considered. Although the 
linear relation between the source strength 
m and the ship surface 


y = tf (x,z) (11) 
is 
ns x 2 (12) 


the actual double model ship body streamline 
should be obtained by solving 


2 =~ (13) 


through the stagnation point where u and v are 
the velocity components of the total velocity 

caused by the double model source distribution 
m and the uniform flow relative to the ship. 


In the linear ship wave theory a smooth 
source distribution produces two systems of 
regular ship waves: the bow and the stern 
waves starting from the bow and stern, 
respectively, represented!! as 


n/2 
S(x,y) = if A(8) exp {i Ss} (x,y,9)} a8 (14) 
-7/2 


aS Seco {(x-x,) cosé + y siné} (15) 


x, the x coordinate of the bow or stern 


g = acceleration due to gravity 


A(8) = amplitude function which is a 
function of source distribution m 


The regular wave ¢ is the solution of 
linear ship wave theory far from the ship. 
Actually, it is easy to see.that the exponential 
function satisfies both the Laplace equation 
and the linear free-surface boundary condition 
for any value of xj. Havelock interpreted 
Equation (14) by discussing the integrand as 
elementary waves,! »13 i.e., the regular waves 
are aggregates of elementary waves starting from 
the bow and stern of a ship. The normal 
direction of each elementary wave crest is 
n= i cosé + j sind. 


Because the local disturbance of the double 
model decays rapidly away from the ship, even 
in nonlinear theory in the far field, the 
regular wave should be of the form of Equation 
(14) with possibly different values of xj and 


A(6). When the integral of ¢ is Oe by 


the method of the stationary phase~’ by taking 
roots of 
3 
su (x+y.0) = 0 (16) 
or 
2 tan*e + = tané +1 =0 (17) 


in general two values of 6 are obtained for a 
given value of each x and y, satisfying, 


rT > 8? (18) 


Furthermore, when 


Xx. g2 
ly 
ar 
Jtano| = 8 2 
or 
6 = 35 deg 


and the vaiue of the integral of Equation (14) 
for any x,y in x/|y [282 can be evaluated by 

the stationary phase method as the sum of two 
waves: transverse waves (0 deg < |6| < 35 deg) 
and divergent waves (35 deg < To| < 90 deg). 


Except for the amplitude function exactly 
the same result for the relation of x,y and 6 
as in the linear theory can be obtained from 
Equation (8) by substituting u = 1 and v = 0. 
Thus, this means that the energy of each 
elementary wave propagates on the uniform flow 
along a straight line ray. The rays are con- 
fined in Sei > 8%, both of the elementary 
waves at 6 = 0 and 6 = 90 deg correspond to the 
ray y /x = 0, and 6 = 35 deg corresponds with 
the ray x/|y| = 82, One ray between these two 
rays corresponds to two elementary waves where 
one is transversel and the other is divergent. 


When the velocity is not uniform due to the 
flow perturbation caused by a ship, the ray is 
not straight but curved near the ship as 
shown in Figure 1. It is known that wave energy 
flux divided by the relative frequency with 
respect to the coordinates for which the fluid 
velocity is zero, is conserved along the 
curved ray cupene If the coordinates are fixed 
in space, then the ship waves are unsteady 
relative to the fixed coordinates and the 
frequency is 


Groiie U k cosé (19) 


-0.1175; 0.08 


O00 = -0.66892 
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Figure 1 - Ray Paths for a Wigley Hull 
(b = 0.1, h = 0.0625) 


because only the transformation of coordinates 
(x,y) to (x - Ut, y) is needed to get the un- 
steady flow. Strictly speaking the relative 
frequency is 


8.02 k u cosé + k v sin6 = k U cos6 
but for our discussion, the approximate value 
serves the purpose. This means that for a 
uniform flow the wave energy is constant 
along the ray because for a uniform flow u = U, 
v=0, k= g/ (U2 cos26), and 8 is constant 
along the straight ray. However, the wave 
energy is not constant along the curved ray in 
nonuniform flow, but is dependent upon the 
local velocity components and 6 because the 
relative frequency is a function of the local 
velocity and 6 as in Equations (5) and (19) and 
@ changes along the curved ray. The wave 
energy will change considerably along the 
curved ray near the bow or stern because u and 
v change to zero at the wave source, the 
stagnation point, while in the far field it 
will be constant along the ray as in linear 
theory. Near the stagnation point, the wave 
number increases according to Equation (5) and 
the wave energy also increases due to the energy 
conservation law.8 Thus, the wave may break near 
the stagnation point. 


Because wave amplitude is so difficult to 
obtain from the ray theory2, the linear theory 
has been considered as an approximation.1 A 
good result was obtained for a two dimensional 
flow problem.9 However, extreme care is needed 
in the three dimensional theory because 8 chang- 
es by a large amount along the ray, near the 
stagnation point, and the wave energy depends 
upon 6. The matching amplitude function of ray 
theory with the values from linear theory should 
be done at the far field where 6 and the direct- 
ion of the ray are, respecitvely, identical for 
hoth cases for each elementary wave. In addi- 
tion, the initial condition has to be taken in 
the neighborhood, but not exactly at, the stag- 
national point because, although the stagnation 
point is the wave source, at the stagnation 
point the ray equations are indeterminant, u 
and v being zero. However, because of the large 
change of u and v near the stagnation point, the 
ray path is very sensitive to the initial values 
of x, y and 6; the nearer to the stagnation 
point, the more sensitive. Therefore, the iden- 
tification of ray paths should be correlated 
with the values of 6 at infinity. Then all the 
ray paths can be properly and uniquely identified 
by 80. 


Since perturbations of the ship decay 
rapidly away from the ship, 6 also rapidly 
approaches 60. The relation between the initial 
value and the value at infinity of 6 has little 
meaning, although it was misunderstood before!» 
because 6 changes very rapidly near the stag- 


nation point. 


Rays of Ship Waves and Ship Boundaries 


Ship waves created by a smooth ship hull 
propagate as regular bow and stern waves from 
the bow and stern stagnation points to infinity 
along rays. Because a ray carries the wave 
energy it cannot penetrate the ship surface. If 
the linear free-surface condition is considered 
with the exact hull boundary condition, as has 


been popular in recent ship wave analysis , 
the straight rays which pass through the hull 
boundary have to be considered. For a ray 
theory with the exact hull boundary condition, 
the ray is not allowed to penetrate the ship. 
In fact, when the initial wave crest touches 
the ship boundary it can be proved that the 
ray of such a wave grazes along the ship 
boundary without penetrating the ship boundary. 


In Equation (8) when the wave crest 
touches the hull 


u cosé + v siné = 0 (20) 
because 6 is the angle between the normal to 
the crest and to the*axis, and the velocity 


normal to the ship hull is zero on the hull. 
Equation (8) then becomes 


yar’ 
dx u 


showing that the ray touches the ship hull 
streamline from Equation (13). 


When the wave crest touches the hull, 
from Equations (11) and (13) 


Vv 
a f (21) 
and from Equations (20) and (21) 


tané = (22) 


we 
x 


Differentiating Equation (22) with respect to 
x along the ship hull 


£ 
+: i (23) 
geal 
x 


dowever, inserting Equation (20) into Equation 
(9) yields 


du , dv 2, du dv 
40 cote ox + a cot 6 ay - coté ay (24) 


u(l - = coté) 


Differentiating Equation (21) with respect to 
x along the hull yields 


v du av v du dy dv oy 


_ uw ox ' ox 7 u dy Ox dy ax (25) 


Vv 
u 12} 


@ 
* 


From Equations (20), (21), (22), and (25) 
it can be shown that Equations (23) and (24) 
are equivalent. That is, when the wave crest 
touches the ship boundary, the ray equations 
and the ship hull streamline equations are 
equivalent. 


When the wave crest is perpendicular to 
the ship hull 


= tané (26) 


El< 


If Equation (26) is inserted into Equation (8) 


dy . ysl 
aie tan 7 (27) 


The ray also touches the ship initially. 
However, Equation (9) is not compatible with 
Equation (26) on the hull. This can be proved 
in a similar way as follows: 


Differentiating Equation (26) with respect 
to x along the ship hull 


£ 
chs pene SE (28) 


dx f 2 cewl 
x 
Inserting Equation (26) into Equation (9) gives 
dé 2 tané (= + tané a) - a ee + tand x) 
— 5 ax dy ay dy 
ahs u (1 + < tané) 


<1 


From Equations (25), (27), and (29), noting 


On oe 
ay ox 
—2£ 
eee (30) 
dx £ ar i 
x 


Equations (28) and (30) are compatible only 
when f,, = 0, or fy = constant. 


That is, only when the ship is a flat plate 
does the ray of the wave, whose crest is perpen- 
dicular to the ship, follow the ship boundary. 
This means that when the ship bow is like a 
wedge, the ray of the bow wave, whose crest is 
perpendicular to the wedge surface, initially 
follows the wedge surface. 


When the ray equations are solved numer- 
ically by the Runge-Kutta method, with initial 
values near the bow stagnation point, the ray 
touches and follows the ship boundary at 


T 
6 So te (31) 


where o is the half entrance angle of the ship 
bow and 6, denotes the initial value of 6. 


When 6; increases from this value the ray 
moves gradually away from the ship as shows in 
Figures 1 through 4. The rays are curved near 
the ship but at far downstream they are straight 
and the ray direction becomes exactly the same 
function of 8 as the linear theory. Thus at 
infinity the rays are inside |dx/dy| = 8%, 
However, when 6. approaches zero, the curved 
ray near the ship gradually approaches the 
ship, and eventually crosses the ship boundary, 
as in Figures 1 through 4. Here the wave 
reflection should be considered to prevent the 
Tay penetration of the ship hull. 
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Figure 4 - Ray Paths for a Wigley Hul] (b = 0.5, h = 0.0625) 


Ray Reflection 


There are no waves coming from places 
other than the ship bow or the ship stern. In 
the ray theory, the flow field perturbed by the 
ship deflects the ray path starting from the 
ship bow toward the ship. Thus, some rays 
of bow waves impinge oa the ship hull. 

Because the ray theory is for small Froude 
humbers and the wave phenomenon is considered 
only on the free surface, it is reasonable to 
consider only reflected waves, as in geometrical 
optics, neglecting transmitted waves when the 
oncoming waves impinge on the ship hull. 


Let the ray at the wave angle 6 on the 
ship boundary (x,y) be reflected to 8; and 
fy = tan®, as in Figure 5. 


Then, 


or 


6 = 26 - 86 (32) 
r ° 


Tigure 5 - Angles of Waves Reflecting on 
the Ship Hull 5 


Whenever the ray intersects the ship 
boundary at (x,y) with angle 8, then the value 
of 6 at (x,y) will be changed to Or = 28, - 8, 
and x,y,8 are continuously calculated by the 
Runge-Kutta's method. Then the ray will reflect 
as in Figures 1 through 4. Here if ® is zero, 
it is easily seen in Equation (8) that the 
impinging ray angle (H) changes to -(H) for the 
reflected ray angle. This fact can be easily 
shown to be true even for 8 # O by just the 
rotation of the coordinate system. 


Numerical Experiments of Ray Paths and 


Free Surface Shock Waves 


Because the ray equations can be solved 
only numerically, careful numerical experiments 
with the ray equations may give valuable in- 
formation. For simplicity, the Wigley hull 
source distribution with 


1 df 
1 era | 
2m dx 
where 
£, = 2nb (1 - 8) (i @) (33) 


is considered for various numbers of b and h 
which are related to the hull beam and draft, 
respectively. The actual hull shape corre- 


sponding to the source Equation (33), is ob- 
tained by plotting the body streamline passing 
through stagnation points as the solution of 


where u and v corresponding to Equation (32), 
are shown in Appendix A together with u.,u, 
and v, for the ray equation. Equations (8)% 
(9), and (10), together with the streamline 
equation, are solved by the Runge-Kutta method 
with initial conditions (x,y,®) near the stag- 
nation points with various values of 8. 


Many ray paths both reflecting and non- 
reflecting from the surfaces of various Wigley 
hulls are shown in Figures 1 through 4. These 
paths were computed by a high speed Burroughs 
computer at David Taylor Naval Ship Research 
and Development Center. The reflection condi- 
tion is incorporated in the high-speed computa- 
tion with a routine to find the intersection of 
the ray and the ship boundary which is pre- 
calculated and saved in the memory. The step 
sizes of integrations and interpolation were 
determined after many numerical tests, and the 
shown results are considered to be reasonably 
accurate. For a given initial condition, the 
solution is stable and converges well. 


In Figures 1 through 4, some common fea- 
tures of rays can be drawn as follows. The 
rays in - 1/2 < 6. < O far behind the ship 
behave like rays of linear theory except wave 
phases are advanced in the ray theory and those 
near’ 6, = O are rays propagating from the ship 
bow and reflecting from those of ship hull. 

The rays near the ship are very different from 
the linear theory as Inui and Kajitani* pointed 
out. The curved rays from the bow have a far 
larger slope than those of linear theory and 
the phase of each ray is considerably advanced. 
The magnitude of the phase difference is more 
sensitive to the beam-length ratio and the 
draft-length ratio than the magnitude of the 
ray slope. 


When 6, = 35 deg, the ray will be the 
outermost ray and ine py angle at ~ will be 
approximately tan” 8-2 as in the linear theory 
and there is the corresponding initial value of 
6 or 64 near the bow, or the origin. However, 
the 6; which is corresponding to a single value 
of 8. is very sensitive to small changes of x 
and y near the origin. At a fixed point near 
the origin there exists a unique correspondence 
between 6, and 6,. 


When from the 64 which is corresponding to 
6, = 35 deg, the initial value of 614 increases, 
8. also increases and the ray angle decreases. 
In general, when 0j = 0, 8. is still a negative 
value. When 6,4 increases further, 68 approaches 
zero and the ray pass is very close to the 
stern. At this point, there exists a certain 
increment of 6; which makes the ray barely 
touch the ship stern, at 04 = 649, When 64 
slightly increases from 84,4» the ray reflects 
from the ship hull near the stern. With the 
increment of 6; the reflection point moves 
toward the bow. When 0 = 84) the ray once 
reflected touches the stern again. When 6; 
increases further from §i], the ray reflects 


twice from the ship hull. In this way, further 
increment of 0; makes the ray reflect from the 
ship hull three, four ... times. However, at 
dj mear the value of a, the ray tries to pene- 
trate the ship hull at the starting point of 
the bow. This cannot be allowed because this 
kind of ray should come from outside of the 
ship. Let the border point of 6; be 6jy- 

This means that rays of initial value of 6 
between 8j, < 989i < Siu reflect from the ship 
hull. As is clear in Figures 1 through 4, 

in general, all the rays before reflection do 
not intersect each other. However, reflected 
rays intersect other rays. The once reflected 
rays intersect not only with each ray once 
reflected from ship boundary points close to 
each other, but also with at least one ray 
before reflection. 


In the stationary phase, each ray has an 
amplitude. Likewise in the ray theory, each 
ray carries it's energy. The reflected ray may 
have approximately the same energy as the ray 
at 6; = 8j0 or 6 5 0 where the amplitude 
function of the linear theory is in general 
more significant than amplitudes of the other 
values of 60. Because the phase must be 
approximately close to each other for the waves 
near 6 = 0, the wave height of the once re- 
flected ray may be close to three times that 
of the transversal wave for 6. = 0. When the 
envelope of the once reflected waves is drawn, 
the domain bounded by the ship surface and 
the envelope, denoted by Dm, must be distinctly 
different from other domains because in Dm 
there are not only once reflected rays but also 
twic or multi-reflected rays on which more than 
three reflected rays intersect by an argument 
similar to that used for the once reflected 
rays. The envelope of the once reflected rays 
behaves like the shock front which was observed 
in Japan. 


In general, a line is called a caustic when 


on one side of the line one can find a continuous 


distribution of rays but not on the other side, 
and along the caustic the wave slope is found 
to be large. The wave near the ray angle 60 = 
35 deg is the caustic, and the wave height near 
the caustic far downstream can be obtained by 
an application of the Airy function in the 
linear theory. The envelope of the once re- 
flected rays may be a kind of caustic formed 

by the refracted bow wave rays due to the non- 
uniform flow perturbed by the ship. Shen, 
Meyer, and Keller? studied such caustics caused 
by the sloping beach of channels and around 
islands. Thus, the additional caustic of ship 
waves may be called the second caustic of ship 
waves, and should not be confused with the first 
caustic which is the known caustic at 6, = 35 
deg. 


Second Caustic of Ship Waves of Various Ships 


For the Wigley hull, several different 
values of parameters b and h were taken to 
find their effect on the second caustic. In 
addition, the effect of a hulbous bow on the 
second caustic was considered. The most 
distinguishable physical characteristic of the 
second caustic is its distance from the ship 
hull. This is related to the distance of 


reflected rays from the ship hull. If the 
number of reflection rays increases, or if 64 
increases from 8ij9, the maximum distance be- 
tween the ray and the ship hull decreases. 

The distance, before or after the ray reflection, 
approximately behaves like a sine curve. The 
maximum distance between the ray before the 
first reflection and the ship hull divided 

by the x coordinate of the point of the first 
reflection a/xy is plotted in Figure 6 for 
various ships. The value of a/x, for different 
values of §{j are approximately the same for a 
given huli and are related to the area between 
the second caustic and the snip hull where there 
may be breaking waves or turbulent waves. Thus, 
if the area is large, viscous dissipation of 
energy becomes large. Accordingly, the measured 
momentum loss behind the ship for the breaking 
waves’ becomes large. 


The values of a/x,, increase with increas- 
ing beam-length ratio. However, the most 
interesting part is the effect of the bulbous 
bow.!0 When the bulb size is increased or the 
doublet strength is increased the curvature of 
the ray near the bow becomes less, although the 
streamline near the bow is such that the en- 
trance angle is slightly large. The values of 
a/x, decrease with increasing bulb size, and 
eventually the ray for 6, = 0 cannot propagate 
without penetrating the ship hull at the begin- 
ning. That is, there is no reflecting ray 
coming out of the bow with a bulb of the proper 
gize, as shown in Figure 7. 
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Figure 6 - Width of the Second Caustic 
a/x, for Wigley Hulls 
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Figure 7 - Rays of the Wigley Ship (b = 0.2, h = 0.0625) with Bulb 
(radius, aS 0.0214h, depth, 2) = 0.5h) 


These phenomena associated with the second 
caustic are exactly the same as those of the 
"free surface shock wave" which was observed 
in the towing tank. 


Finally it should be pointed out that the 
ray of the wave whose crest is perpendicular 
to the hull surface follows exactly along the 
flat surface as was proven by Equations (28) 
and (30) and the following paragraphs. This 
means that rays originating from the vertex 
of a wedge very likely never reflect on the 
wedge surface. This is because rays near the 
ray which follows the ship surface did not 
intersect each other near 6; = -1/2 + a, or 
®. = - 1/2. Rays near 6, = 0 impinge on the 
ship surface due to the effect of the water- 
line curvature of the ship. Therefore, if a 
ship has a wedge bow, the second caustic must 
be found near or behind the shoulder. Be- 
cause the first caustic near the bow may be 
very prominent and both waves near the first 
caustic and waves near the stagnation point 
break, careful experimental analysis and more 
theoretical study of the bow near field may be 
needed. 


Wave Amplitude and Phase 


Because the perturbation due to a ship, 
or both regular waves and local disturbances, 
decays at far field, the linear theory must 
hold in the far field. In particular, the 
wave resistance can be calculated from the 
energy passing through a vertical plane x= 
constant far downstream; the linear theory 
which is properly matched to the near field of 
the ship will be used for calculating the wave 
resistance. As explained in Equation (14) and 
the following paragraphs, the expression for 
regular waves far downstream is known to be of 
the form of Equation (14) where the amplitude 
function 


A(8) = P(8) + 4 Q(8) (34) 


may be taken as an approximation from the linear 
theory but the phase difference xj must be ob- 
tained from matching with the near field. If 
the near field is also expressed by the linear 
theory 

1 


P+iQ=T ky i i 


-h o (35) 
2 
Les fan ik, x secé 
dxdz fe sec 0 e 


where f, is the derivative of Equation (11) with 
respect to x. When the inner integrand of 
Equation (35) is integrated with respect to x, 
the value with the limit x = 1 will form stern 
waves and the value with the limit x = 0 will 
form bow waves.!1 Then the bow waves can be 
represented by 


m/2 
Go f A, (8) exp {i s, (0)} de (36) 
-1/2 
where 
A, (0) =P, +4, (37) 


When the phase s is computed from Equation 
(10) along with the ray path from Equations 
(8) and (9) considering that s = 0 at the bow 
near the origin, there are two results different 
from those of linear theory: (1) the ray path 
is deflected as if the elementary wave of the 
linear theory started from (x;,0) not from the 
origin, and (2) the phase change denoted by 
As should be considered. That is, the 


equivalent linear elementary wave may be 
written as 


2 
(6) exp {i kmersecso 
Ay ° 1 


{(x-x,) cosé + y siné} ar pt ko As] 


where x; is obtained as an intersection of the 
tangent to the ray at ~ and the x axis and 


2 
ko ds = kos - ko sec 6 {(x-x,) cosé + y sino} 
(39) 
The value of 
S, (6) = As - yy secé (40) 


can be obtained at any point along the ray. 

In general, x, is negative and s?2 (8) is 
positive meaning that the bow wave phase in the 
ray theory is larger or more advanced than the 
phase of the linear theory. This fact has long 
been observed in experiments in towing tanks. 


The advancement of wave phase is computed 
for various ships and the values of s?2 and x] 
at x = 2 are shown in Figures 8 through 10. 
When the beam-length ratio increases and/or 
the draft-length ratio increases, the values 
of s? increase and the values of x, decrease 
for all values of 6,. As compared with the 
increment of the slopes of rays near the ship, 
the increment of the phase angle is more 
sensitive to the beam- and/or draft-to-length 
ratios. 


The most interesting phenomenon about the 
phase difference is in regard to the bulbous 
bow.!0 That is, the phase differences for hulls 
with and without bulbs are almost the same even 
with a considerably larger bulb. In the past, 
because of the observed phase difference of 
ship waves, the bulb was located far forward 
to obtain good bow wave cancellation.! 
According to the present analysis, if there is 
no other reason, the bulb position need not be 
far forward. Because the nonuniform flow cre- 
ated by the ship is much more significant than 
that of the bulb, as far as phase change is 
concerned, both the ship bow waves (in general, 
positive sine waves) and the bulb waves 
(negative sine waves) propagate through the 
same region and cancel each other. 


As for the amplitude function, although 
it was shown by Doctors and Dagan? that ray 
theory produced the best result for a two- 
dimensional submerged body even though they 
used a linear amplitude function, the surface 
piercing three-dimensional case may be quite 
different. The amplitude function is mainly 
related to the singularity strength which sat- 
isfies the ship hull boundary condition and 
some improvement might result by considering the 
sheltering effect. However, in the present 
study, the linear amplitude function is used to 
simplify the problem, showing the effect of 
curved rays. 


10 


Wave Resistance 
Ye 


If all the elementary waves are assumed to 
be propagated without reflection, the amplitude 
function and the phase difference studied in 
the previous sections will supply enough infor- 
mation for the calculation of wave resistance. 
Because at far downstream, the wave height 
may be considered linear, the Havelock wave 
resistance formula!> may be used for the waves 
represented by Equations (36) through (40), 
considering that the wave with the changed 


phase Sop (8) has the amplitude 
ikos2, (6) 
(Py +i a e 
w/2 
R 
=C = (41) 
22. w 
2/2 UL S19 
ik s (6) 
(P, +iQ)e as cos26d8 
b b 
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ae 16m ko ej 
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w = width of the towing tank nondimension- 
alized by ship length L 
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Figure 8 - The y Coordinate of Ray at x = 2, Phase Function S7(8) and the Starting Point of 
Elementary Wave, x) of a Wigley Hull (b = 0.2) 
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Figure 10 - S,(8.)> and -x 


the phase difference of each elementary wave 
does not change bow wave resistance from the 
linear theory. 


The stern wave amplitude is exactly the 
same as the linear value 


As (8) = Ee ap iQ. 


but the phase difference s2, (8) should be com- 
puted by the ray theory together with the ray 
path which is shown in Figure 1l. Then it is 
also obvious that the stern wave resistance is 
the same as the linear stern wave resistance. 
The total wave resistance may be obtained simi- 
larly by considering that the total wave ampli- 
tude is 


iks,, (6) 
(P, + iQ)e ap Ge. SP fi) _)) 
b b s s (44) 
se (So¢ (8) - secé) 


Here, the bow and stern wave interaction ap- 
pears in the wave resistance. That is, only 
the interaction term changes due to the phase 
change caused by the nonuniform flow. This 
fact is exactly the same as in two-dimensional 
theory. 


The actual computation of wave resistance 
is performed by the Srettensky formula using 
the relation 


1 2 (45) 


tan. 


and the corresponding values of s7 are obtained 
by interpolation. When a portion of elementary 
waves near $8. = 0 is reflected from the ship 
hull, the larger part of the energy in this 
portion of elementary waves will be dissipated 
by breaking waves, and the wave resistance will 
decrease, but the momentum loss due to breaking 
waves will increase. If such energy was con- 
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of a Wigley Hull (b = 0.2, h = 0.0625) with and without 
bulb (x, = 0.0285h, 2515 0.7h) 


considered to be totally missing in the wave 
resistance, j in the Srettensky formula for bow 
waves would start not from zero but from a 
certain number minimum J = J such that 


where 6.) is the value of 8. from which bow 


waves start to reflect. If 8.0] = Oa5 Ga = Ik, 
and €. = 2 for J > Jj, and if 6; > 6,,~! 
GAS for J > Jj. me 


The result is shown in Figure 12 where a 
considerable shift of the phase of hump and 
hollow of the wave resistance due to bow 
sterm wave interaction is noticeable. 
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Figure 11 - Ray Paths of Stern Waves of 
a Wigley Hull (b = 0.2, h = 0.0625) 
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Figure 12 - Wave Resistance of the 
Wigley Ship (b = 0.2, h = 0.0625) 


Discussion 


Recently Rasece.? obtained a dispersion 
relation from a low speed free-surface boun- 
dary condition which was a slightly modified 
version of the one Babal8 used. Eggers showed 
that there was a small region near the stag- 
nation point where the wave number became 
negative; since this was not permissible fora 
wave and it could be interpreted to mean that 
there was no wave in the region. He suggested 
that his form might help alleviate the sen- 
sitivity of the initial condition to the ray 
paths. According to the energy conservation 
law8 of the wave propagating through non- 
uniform flow the wave energy flux is proportion- 
al to the wave number along the ray tube. 

When the wave number along a ray is considered, 
although the Keller dispersion relation has an 
infinite wave number only at the stagnation 
point, the Eggers dispersion relation has an 
infinite wave number in the flow near the bow. 
Thus, near such a singular point or a singular 
line, waves may break and the present theory 
cannot be applied. When Eggers equation was 
incorporated into the ray equation in the 
present ray computer program it was found that 
the ray path was still very sensitive to a 
small change of initial value (x,y,8) and the 
curved ray still reflected from the ship hull. 


In the experiments’ conducted at Tokyo 
University, the second caustic can be noticed, 
and near the second caustic the flow field is 
violently different from the linear theory. 
However, it seems that extreme caution is 
needed to distinguish the second caustic from 
the first caustic. The case of experiments 
with wedges is interesting because, accoraing 
to the present theory, there cannot be a second 
caustic on the wedge although the first caustic 
should be there. However, the wave phases of 
the wedge should be advanced and ‘quite sensitive 
to the draft- and beam-length ratios due to the 
nonuniform flow caused by the wedge. On the 
other hand, the present theory is still not 
exact although it takes into account the effect 
on the propagation of water waves of nonuniform 


3) 


flow due to the double model. Thus, the non- 
linear effect of water waves is not totally 
analyzed here. Nevertheless, the present 
theory supplies a great deal of hope to an 
entirely different approach to the ship wave 
theory - the ray theory. 


Although the strictly linear amplitude 
function is used for simplicity in the present 
study, a slightly improved amplitude function 
may be easily incorporated by adding the 
sheltering effect or other effects. However, 
it should be noted here that, in any case, the 
bow wave resistance and the stern wave resis- 
tance are the same as the results without the 
ray theory and the effect of the ray theory 
would appear as a shift of the phase of hump 
and hollow in the total wave resistance. There- 
fore, the computational results would not match 
the results of towing tank experiments unless 
the viscous boundary layer effect on the, gtern 
wave or some other effect is considered. 
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Appendix A 
For the computation of ray paths of a 
ship, the flow velocity and its derivatives on 


(x,y,0), U,V,uy, V,uy are needed. A Wigley 
hull has the doublé model source distribution 
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